Creep via dynamical functional renormalization group 
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We study a D-dimensional interface driven in a disordered medium. We derive finite temperature 
and velocity functional renormalization group (FRG) equations, valid in a e = 4 — D expansion. 
These equations allow in principle for a complete study of the the velocity versus applied force 
characteristics. We focus here on the creep regime at finite temperature and small velocity. We show 
how our FRG approach gives the form of the v — f characteristics in this regime, and in particular 
the creep exponent, obtained previously only through phenomenological scaling arguments. 



Lines and surfaces can exhibit a remarkable variety of 
complex phenomena when they are driven by an external 
force over a disordered substrate. This situation is en- 
countered in numerous physical systems such as growth 
phenomena 0^], magnetic domain wall motion ||, fluid 
invasion of porous media jj] , vortex systems [0 , charge 
density waves (CDW) |(| or Wigner crystals 0. 

In all these cases there is a competition between the 
elastic energy, which tends to keep the system straight, 
and the disorder. In addition to making the static inter- 
face rough, the disorder has also important consequences 
for the dynamics. At zero temperature T = 0, the inter- 
face remains pinned until a critical force f c is reached, 
whereas for large drive, it moves with a velocity propor- 
tional to the force, the disorder being averaged by the 
fast motion. Obtaining theoretically the velocity v ver- 
sus applied force / characteristics is a very challenging 
problem since such a v — f curve is directly measurable, 
and is in most cases one of the most important physical 
properties (e.g. the transport properties in vortex, CDW 
or Wigner crystal systems). 

A very fruitful approach is to cast the depinning tran- 
sition of an elastic system in the general framework of 
critical phenomena, with the velocity as an order param- 
eter A functional renormalization group (FRG) 
treatment of the equation of motion ^,|) allows to obtain 
the various critical exponents characterizing the depin- 
ning transition. 

At finite temperature, the situation is more complex 
since the interface can move by thermal activation below 
the zero temperature threshold force. This leads to a 
rounding of the depinning transition (see e.g. | fl2|| ) when 
/ ~ f c , and to a creep motion of the interface when / <C 
f c . Indeed it was realized that, due to the glassy nature 
of the static pinned interface, a moving interface would 
have to overcome divergent barriers as the applied force 
is reduced [|l3|— p~7[| in contrast to much older theories that 
assumed that motion occured only through finite barri- 
ers Jl8| . This lead to the proposal of a phenomenological 
theory of creep resulting in a highly non- linear v — f char- 
acteristics of the form v cx exp {—{U c /T)(f c /f )^). Using 
scaling arguments, and assuming that relevant barriers 
for the dynamics scale with the same exponent as the 
energies of metastable static configurations, one obtains 



M = {D — 2 + 2C)/(2 - C) where D is the dimension 
of the interface and £ the static roughening exponent of 
the interface. This remarkable formula relates dynamical 
properties to purely static quantities. The creep theory 
has been extremely successful in explaining various phys- 
ical phenomena, in particular vortex systems Q , and the 
quantitative relations between the exponents has recently 
been verified by experiments on magnetic domain walls 




FIG. 1. Typical force- velocity curves, 
with relative transverse displacement u. 



Inset: an interface 



Unfortunately, up to now, the creep formula has re- 
mained phenomenological ||l9|| and relies on a certain 
number of mostly unverified assumptions |20[ . On the 
other hand, extending to finite temperatures and veloci- 
ties the FRG techniques used at T = has proved very 
challenging, although some interesting results could be 
obtained for the periodic systems driven within their in- 
ternal space pH-pit . We address this issue in the present 
paper. We give a generalization of the FRG equations 
valid for arbitrary temperature and velocity. These equa- 
tions allow in principle for a complete study of the ther- 
mal effects on the v — f characteristics. We show in this 
letter how they can be used to obtain the creep behav- 
ior and to derive the creep exponent from first principle 
(i.e. the equation of motion). We compare our findings 
with the phenomenological theory. A more refined study 
of the flow, a precise calculation of prefactors, and other 
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consequences will be presented in a longer publication 
II- 

We consider a D-dimensional manifold without over- 
hangs described by a height function u, as shown in Fig- 
ure , embedded in a space of dimension d. If D = d this 
also describes periodic systems The system obeys 

the equation of motion 

{T]d t ~ cV 2 )u rt =f-rjv + F(r, vt + u rt ) + ( rt (1) 

where rj is the friction coefficient and c is the elastic con- 
stant. For simplicity we assume here isotropic elasticity, 
but generalization to less symmetric elastic tensors can 
easily be done. Thermal fluctuations are described by 
the Langevin force (CrtCr'f) = 2r)T5 D {r - r')6(t - t') 
(the brackets denote thermal averages). Disorder gives 
rise to a random force characterized by F(r, u)F(r' ,u') = 
5 D (r — r')A(u — u') (the overbar denotes disorder aver- 
ages). 
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FIG. 2. Typical correlators for different cases of disorder 
and the corresponding static roughening exponent 

Depending on the precise microscopic nature of the 
problem, the correlator A can have different forms ||] as 
shown on Figure . The random field case occurs when A 
is short range, the disorder is random bond when the force 
is the derivative of a potential which has short range cor- 
relations, and finally, when the system is periodic, such 
as vortex lattices, A is periodic [^. The total height of 
the interface is vt + u rt , and the condition (d t u r t) = 
fixes self-consistently a relation between the average ve- 
locity v and the force /. In the absence of disorder, we 
are left with the trivial friction law j = r\v. In the fol- 
lowing we denote by / = / — r\v the deviation from this 
law. 

To make a systematic perturbation expansion and 
renormalization out of this stochastic equation, we fol- 
low the standard method of dynamic field theory [ p7|j28| 
which yields an action S(u,u), where u is the response 



field. Within this formalism, a perturbative expansion 
in disorder A allows to generalize |^3| the large veloc- 
ity expansion |29|]3C| ] of the characteristics f(v) to finite 
temperatures. But the main interest of this formalism is 
to allow to go beyond perturbation theory and to build 
the FRG equations at finite temperatures and velocities. 
The procedure is similar to the one used in p3| ], so we 
just give an outline here and leave technical details for 
p5j. We expand the weight e~ s up to second order in 
A, knowing that / = 0(A). We then integrate over the 
"fast" modes of the fields u and u having Fourier compo- 
nents Ae~ l < q < A, where 1/A is a short distance cut- 
off. Reexponentiation gives an effective action similar to 
the original one, but with renormalized parameters, once 
we restore the original cutoff by rescaling the fields as 
u'(r,t) — e~^ l u(e l r,e zl t) and u'(r,t) = e~ al u(e L r, e zl t). 
The elasticity c is not renormalized, for statistical tilt 
symmetry reasons . This fixes a such that the triv- 
ial scaling of c is zero. We are then left with two free 
parameters (which can be chosen /-dependent if needed): 
C, the roughening exponent of the interface relating longi- 
tudinal and transverse lengths u ~ r^, and z, the dynam- 
ical exponent relating space and time t ~ r z . Denoting 
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A by A, 



by T, Jgf by A, keeping A fixed 

(Sd is the surface of the sphere divided by (2ir) D ) and 
performing the rescalings, we obtain the following flow 
equations in a e — 4 — D expansion 

9 ; A(m) = (e - 20A(u) + C«A'(«) + TA"(u) 

e - " [A'» (A((V - s)A) - A(u + (s' - s)A)) 

s>0,s'>0 

-A'(u-s'X)A'(u + sX) 

+A'((s' + s)A) (A> - s'X) - A'(u + sX))] (2) 



rjv 



di\n\ = 2-(- [ e - s sA"{sX) 

Js>0 

cUnT = e-2-2C+ f e - s sXA'"(sX) 
Js>a 

dJ^(2-0f + cA 2 [ e~ s A'(sX) 

Js>0 

dic = 
di In v = z — C 



(3) 
(4) 

(5) 

(6) 
(7) 



These equations ]3l{ | allow in principle for a complete de- 
scription of the properties of the interface at finite tem- 
peratures and velocities. 

For the deterministic (T = 0) static (v = 0) fixed point, 
the flow equations are greatly simplified since T\ = 
and Xi = 0. One recovers the usual formulas derived 
in the context of the depinning transition via replicas 
HHI or FRG 10. The correlator A/ flows to a fixed 
point A* which depends on the type of disorder. Ex- 
act solutions are known for random field |B[ and peri- 
odic cases [p6|, and asymptotics for random bond. In 



2 



all these cases, the fixed point function has a cusp (A*' 
is discontinuous at the origin), which appears at a finite 
scale l c — Mog(l + _ 3 a"(o) ) an< ^ which is responsible 
for the existence of the pinning force |3 26 1. The scale 
at which the divergence of A"(0) occurs is the Larkin 
length, typical size of a segment of the interface wan- 
dering over the correlation length of the disorder. Its 
expression e lc /A coincides with its standard Q expres- 

2 2 

sion, R c ~ (e-^fjy) 1 ^! when restoring the original A, and 
denoting its range by 77. In the following, we denote by a 

the length -A*'(0+), where @ a ~ 0.4e 2 / 3 (/A ;=0 ) 1/3 
for the random field case and a — ea/6 for the periodic 
case where a is the lattice spacing. 

The flow equations (§-0) at finite but small velocity 
and zero temperature provide a natural derivation of the 
procedure used in || for calculating the critical force and 
the depinning exponents. Since up to the Larkin scale, 
disorder and displacements are small, the flow is trivial. 
After l c , A has reached its fixed point A*. Provided that 
the velocity is small (A/ is smaller than the range of 
A*), our flow equations coincide with those used in || 
since we can replace the r.h.s. of the flow equations by 
the values at 0+ of A*. The scale used to cut the flow 
in P| appears here naturally when the length A; is of 



the order of the range 



As in H, we denote by ly 



the corresponding value of I. Note that can be quite 
different from the original range rj of the disorder. At 
larger scales, one crosses over for I ~ ly to a regime where 
the disorder acts as Langevin forces, and we are left with 
an effective Edwards Wilkinson equation of motion. We 
thus recover [f25f 



fc 



cA 2 A*'(0+) 
' 2-C ' 



t (C-2)/o 



(8) 



This expression gives back, up to an e coefficient, the 
dimensional estimate f c ~ ^5-, since R c is e lc /A and 
ae C'c ^ o.2ery in the random field case and ae^ ic = ea/6 
in the periodic case. 

Let us now analyse the finite temperature case. A finite 
temperature flows to zero with the free energy exponent 
e - 2 - 2C- But the presence of the TA"(u) term in (g) 
is enough to smooth the behavior of the correlator near 
the origin, and it removes the divergence of A"(0). In- 
deed one expects A; to remain analytic at all finite scales 
even if its limit A* has a cusp. More precisely, during 
the flow, — A"(0) increases and one may assume |3j| that 
— Ti A"(0) converges, since the flow equation at the origin 
is^A/(0) = (e-2C)A / (0) + T / A; / (0). One can check that 
its limit is a 2 (with a defined above). Using this prop- 
erty in the T > 0, v — FRG equations for A and T we 
obtain the precise way ^{| the successive derivatives of 
A grow with I, and get 



1(^(0) -aA) 
Ti \ a 



The rounding close to the origin appears at a scale of 
the order of the temperature TJ. The form of A; in this 
regime is given by A* for u = ^ and by the round- 
ing (||) for u <C as shown on Figure . 









r* 







1 



(9) 



FIG. 3. Generation of the cusp: two lengthscales appear, 
the width £j of the rounding, of order Ti, and the width r*j of 
the fixed point function A*. 



This property of the finite temperature equations en- 
ables us to investigate the creep regime. In addition to 
the Larkin scale l c , two obvious lengthscales appear in 
the equations: the scale It for which A; = £j, the width 
of the rounding of A;, and as before, ly for which A; is 
of the order of the range of A*, and above which dis- 
order is washed out. Provided / <C f c , the velocity is 
arbitrarily small, and so is the initial value of A. One 
has thus l c < It < W ■ The equation for the tempera- 
ture introduces in fact an additional lengthscale that we 
discuss later. Ignoring for the time being this additional 
lengthscale allows to distinguish two main regimes above 
the Larkin scale l c . In the thermal regime l c < I < It, 
we can replace the set of equations (||j3) by static ones 
obtained by replacing the r.h.s. occurences of A; by zero, 
and not by 0+ as for the depinning case. The quan- 
tity — T;A"(0) becomes constant until A; ~ In this 
regime A grows rapidly. Physically this corresponds to 
a regime where barriers grow and the system can only 
move through thermal activation over those barriers. 

Beyond It one enters the depinning regime. The dis- 
order is close to its fixed point and the other quantities 
now flow with the depinning flow equations, since the in- 
tegrals over s do no more see the behaviour at the origin 
but the completely different values of A* at 0+. Roughly 
speaking we can in this regime ignore the temperature, 
as far as motion is concerned. Physically this is because 
above the length e' T /A - i? c (/ c //) 1 /( 2 ~ c ), that corre- 
sponds [^H to the criterion A; ~ the external force is 
dominant over the pinning force, and these lengthscales 
do not need thermal activation to move. For I 3> It, 
the temperature is not important for the motion. As for 
T = 0, we cut the flow at the scale ly where a perturba- 
tive expansion in disorder is well defined. We can use it 
to estimate the value of fi v . Since the renormalized tem- 
perature is very small it is enough to do this calculation 
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as if Ti v = 0. Then it becomes obvious that fi v does not 
depend on the initial velocity, force or temperature. 

Solving the flow (|2|jj]) for the three regimes yields a 
relation between initial quantities. Using (||) and replac- 
ing the initial value of / by / since in the creep regime 
/ » rjv, we get 



(2 - C)Sd 



Id 

( e / c /A)^-2( Q , e ic)2 c ^f c 



1+M 



exp 



i 



(10) 



where fj, = (2 
U c gives 



e + 2Q/(2 — C) and a rough estimate of 



c(ae^) 2 (e l <=/A) 2 



{2-e + 2C,)S D (2-e + 2()S L 



(11) 



where C is a constant. We thus recover precisely the 
creep exponent /i, obtained previously only via the phe- 
nomenological barrier arguments. Our derivation allows 
to obtain this value directly from the equation of mo- 
tion without having to do any supplementary assump- 
tion. The fact that we do recover the simple scaling ex- 
ponent /z, is thus also a strong indication that the barriers 
important for the motion do scale with the free energy 
exponent (up to subleading corrections). Although our 
calculation is not expected to give, at that level, a pre- 
cise estimate of the prefactors in the barriers U c it is 
contorting to also recover an expression similar to the 
one derived by scaling arguments, U c ~ cr 2 R^~ 2 , i.e. 
the barriers at the lengthscale of the Larkin length. As 
we showed above the FRG also identifies correctly the 
relevant lengthscales appearing in the simple scaling ap- 
proach of the creep. This simple application of the FRG 
equations then gives a consistent picture for the creep. 

The main interest of the present method is of course 
to allow in principle to go beyond the simple scaling. 
In this respect various questions remain. In particular 
the FRG equation for the temperature leads to an addi- 
tional lengthscale at A/ ~ T^ 2 . Above this lengthscale 
the disorder starts renormalizing the temperauture up- 
wards. Although this effect has no impact on the creep 
exponent pq] it can obviously lead to a modification of 
the barriers (or equivalently of the 1/T term in (|To|)). 
A more precise calculation of this amplitude, taking into 
account the full crossover around I ~ It will be presented 
in |2f| . Given the success of the finite temperature FRG 
for the creep regime, it is also tempting to apply these 
methods to the study of other regimes of the motion. In 
particular, close to the threshold at finite temperature, a 
scaling approach describing the motion as a succession of 
avalanches, and a complementary FRG approach based 
on the flow equations derived here, should help fixing the 
form of the force-velocity curve. 
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